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^P^; ABSTRACT 

I ■ We present high spatial resolution MERLIN 1.4GHz radio observations of two high 

P ; redshift (z-2) sources, RGJ123623 (HDF147) and RGJ123617 (HDF130), selected 

as the brightest radio sources from a sample of submillimetre-faint radio galaxies. 
^ ■ They have starburst classifications from their rest-frame UV spectra. However, their 

radio morphologies are remarkably compact (<80mas and <65mas respectively), de- 
manding that the radio luminosity be dominated by Active Galactic Nuclei (AGN) 
rather than starbursts. Near-IR imaging {HST NICMOS F160W) shows large scale 
^ \ sizes (Ri/2 ~0.75'', diameters ~12kpc) and SED fitting to photometric points (optical 

through the mid-IR) reveals massive (^SxIO^^Mq), old (a few Gyr) stellar popula- 
tions. Both sources have low fiux densities at observed 24/imand are undetected in 
observed 70/im and 850/im, suggesting a low mass of interstellar dust. They are also 
formally undetected in the ultra-deep 2Ms Chandra data, suggesting that any AGN 
activity is likely intrinsically weak. We suggest both galaxies have evolved stellar pop- 
ulations, low star formation rates, and low accretion rates onto massive black holes 



oo 

(N 

in 



(N 

o 

I (IO^'^Mq) whose radio luminosity is weakly beamed (by factors of a few). A cluster- 

like environment has been identified near HDF130 by an over-density of galaxies at 
z = 1.99, reinforcing the claim that clusters lead to more rapid evolution in galaxy 
populations. These observations suggest that high-resolution radio (MERLIN) can be 
a superb diagnostic tool of AGN in the diverse galaxy populations at z~2. 
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1 INTRODUCTION 

The development of deep, mult i wavelength imaging sur- 
veys has shed light on the rapid construction of galaxies 
and the build-up of stellar mass at high redshift (z~2- 
3). The population of galaxies we see at this epoch is 
diverse— a testament to their rapid and significant growth. 
The cosmic star formation density peaks at this redshift 
(z^2) as does the q uasar space density (Fan et al. I I2OOII : 
iBlain et al]l2002l : iRichards et al. 2006). Galaxies at z^2 dis- 
play a wide range in properties, from negligible nuclear 
activity to the most powerful AGN, and with star for- 
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mation rates (SFRs) over 1000 MfT) yr~^ in Submillime- 
tre Galaxies (SMGs; e.g. Blain et all l2002h to as low as 
1 yr~^ (e.^. quiesc ent galaxies; IPaddi et al.l l2005l , 

I2OO4I : ICimatti et al.|[20o3 ). Having also come to light in the 
past decade is the importance of AGN in regulating bulge 
growth fe.g. IPi Ma tteo et al. 2005). This is apparent in the 
local M^-Mbh relation where black hole mass is roughly 
^0.1% of the spheroid's stella r mass (|Magorrian et al.lll998l : 
iKormendv Richstonel Il995l : iGebhardt et al.1 12OOOI ). This 
relationship clearly suggests that black hole mass and stel- 
lar mass build up symbiotically, with a broad range of 
galaxies lying withi n a factor of a f ew of the lo- 
cal M^-M^jf ratio (e . g.. IPeng et al.1 I2OO6I : iMcLure etHI 
I2OO6I : [Alexander et al.1 l2008bi r Moderate-luminosity AGN 



with apparently low accretion efficiencies, ultra-luminous 
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AGN with powerful feedback winds, and heavily obscured, 
Compton thick AGN whose strength is shielded by dust 
in the host galaxy have been identifi ed in galaxies with 



Table 1. Physical Properties of HDF147 and HDF130 



^ at high-z (Brandt fc HasingerllioOsI: 
20071 : IComastril 



Alexander et al.ri2005l . l2008al : IPaddi et al.1 
20041 1. Only through the most systematic searches of deep 



survey data may the full range of galaxy properties be 
clearly understood. 

In an effort to study an active subset of galax- 
ies, our collaboration has focused on the /iJy radio galaxies, 
which are roughly 1/3 AGN, 1/3 submillimetre-bright star 
formers (SMGs) and 1/3 submillimetre- faint star formers 
(submillimetre- faint star forming radio galaxies, SFRGs). 
While most SMGs and SFRGs may have some contribu- 
tion from low- luminosity AGN, their bolometric luminos- 
ity is dominated by dust-generated emission in the far-IR. 
The AGN sources have bolometric luminosities dominated 
by non-thermal radiation from AGN. SFRGs are likely to 
be the hotter dust c ounterparts of the SMG population 
([chapman et al.ll2004l ). the difference in their submillimetre 
properties caused by a shift in the blackbody dust emission 
towards higher temperatures (T^ >50K). There is no other 
obvious physical difference between SFRGs and SMGs, both 
contributing significantly to the z^2 star formation history 
(Casey et al. 2008, in preparation). Here we present de- 
tailed mult i wavelength studies of two galaxies, RG J123623 
(HDF147) and RGJ123617 (HDF130), originally identified 
as SFRGs but postulated here as being radio-bright but in- 
trinsically low- luminosity AGN. 

In ^ of this paper we describe the observations and 
the parent sample from which the two galaxies are drawn. 
The observations include the MERLIN radio survey, VLA 
coverage of HDF, coverage by Spitzer 24/im, 70/im, HST 
photometry and morphology, and X-Ray from CDF-N. In 33] 
we interpret these observations and derive an SED model, 
while in 311 we discuss the selection population and the pos- 
sible physical scenarios that could describe such a set of 
mult i wavelength observations, from a massive yet poorly ac- 
creting black hole whose radio emission is weakly beamed, 
to strongly beamed low- luminosity AGN, and to heav- 
ily obscured Compton thick AGN. Throughout we assume 
Ho=70kms"\ ^^0=0.3 and ^^a=0.7. 



2 SAMPLE AND OBSERVATIONS 

HDF147 and HDF130 (see Table 1 for observation details) 
originate fro m a sample of 18 SF RGs (also referred to 
as 'OFRGs'; 'Chap man et al.l [20041 ). The sample was de- 
fined by all radio-detected galaxies in seven distinct fields 
that were previously targeted in the sub-mm by SCUBA 
([Holland et al.l Il999[ ) but lacked an opti cal counterpart 
brighter than R — 23.5. As described in IChapman et al.l 
([20041 1. the requirement of an optically faint source elimi- 
nates low redshift, moderate- luminosity sources as well as 
high redshift, bright AGN and leaves a proposed population 
of predominantly hot dusty ultra-luminous galaxies at z ^ 
2. The set of SFRGs have much lower submillimetre flux 
dens ities (^0.5 mJy) than the nominal ^-^6.5 mJy of SMGs 
(e.g. IChapman et al.1 12005 ). 

High-resolution observations from the Multi-Element 
Radio Linked Interferometer Network (MERLIN; 
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Z mag 
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24.7±0.5 


H mag 
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0.40 


0.03 
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0.75±0.06 
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Notes — All magnitudes are in the AB-magnitude system. X- 
Ray luminosity is given for the rest frame 0.5-8.0keV band. The 
variable R1/2 denotes the half light radius of the galaxy in H band 
and MERLIN radio imaging respectively. The Sersic Exponent, 
n, refers to the Galfit profile fitting to H-band. The derived 
quantities come from calculations described in JS] 

[Thomassoij Il986h were obtained for these sources as 



described in Muxlow et al.l ([20051 ). One motivation for 



the high angular resolution observations is distinguishing 
the compact emission of an AGN from more diffuse star 
formation, sep arating star forming galaxies from AGN. 
IMuxIow et "all (^2005) discusses the 0.3-0.5'' beamsize MER- 
LIN observations in depth. The rms noise from MERLIN 
data was 5.9 aJy beam~^. The im ages from the VLA 
([Richards et al ] I1999I : iRichardsl liooO ) and MERLIN were 
combined in the image plane and used to make a sensitive 
1.4- GHz map with high positional accuracy (tens of mas) 
allowing high resolution imaging and identification of /xJy 
sources. 

Spitzer observations were made in IRAC 3.6/im, 4. 5Mni, 
5.8/im and 8.0/im and MIPS 70/im (Frave r et aLll2006h and 
24/im bands in the HDF region (PI M. Dickinson; R. Chary 
et al., in prep). This component of observations is used to 
study the galaxies' stellar masses. Optical Imaging of the 
area is from the Hubble Deep Fielctl, in B (F435W), V 



^ Based on observations made with the NASA/ESA Hubble 
Space Telescope, and obtained from the Hubble Legacy Archive, 
which is a collaboration between the Space Telescope Science In- 
stitute (STScI/NASA), the Space Telescope European Coordinat- 
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(F606W), i (F814W), and z (F850LP) bands all in the AB 
magnitude system. As part of the extensive 45 arcmin^ near 
infrared survey coinciding with the HDF and GOODS-N, 
we have obtained NICMOS imaging using NIC3 three orbit 
pointing in F160W of both HDF147 and HDF130 (Conselice 
et aL, in prep). 

The parent sample of SF RGs was also observed spec- 
troscopically with Keck/LRIS (Stei del et al.l l2004': O ke et al.l 
[199 5 ) from 2002-2004. Full details of these rest-U V spectral 
observations are found in iChapman et aP (|2004l ) . In short, 
the galaxies' redshifts were identified by detection of several 
emission features (Lya and CIV in the case of HDF147 and 
Lya and Hell in the case of HDF130). Although their spec- 
tra are faint (see Figure 2 of Chapman et al.), both resemble 
starburst galaxies, so we grouped them originally in the par- 
ent SFRG sample. We attempted to measure individual line 
fluxes from several features that would more clearly classify 
HDF147/HDF130 as starbursts or AGN, but the spectra 
are to o faint to go beyond the analysis of IChapman et ahl 
(I2OO4I ). 

Observations in the X-Ray are from the 2 Ms Chandra 
Deep Field North (CDF-N; Alexander et al. 2003), centered 
on the HDF, in the 0.5-8.0 keV band which achieves an 
aim-point sensitivity of ^1.2 x 10~^^ ergs cm~^ s~^. X- 
Ray obse rvational details of a si milar SMG population are 
found in [Alexander et al ] (2005). At the redshifts of these 
galaxies (~1.96), the observed 0.5-8.0 keV band corresponds 
to res t frame energies 1 .5 - 2 4keV. Neither galaxy is listed 
in the l Alexander et al.1 (|2003h catalog; however, HDF147 is 
marginally detected at a low significance in the narrow 1- 
2keV band (a luminosity of 1.67x10^^ erg s~^ in the ob- 
served 0.5-8 keV band). Assuming an average X-ray spec- 
tral slope r = 1.8, we find the luminosity of HDF147 in 
the rest 0.5-8.0 keV band is Lx(hdfi47) = 4.1x10^^ erg s~-^ 
and the the 3-a upper limit for HDF130 is Lx(i^z:»Fi30) (2- 
10keV)<8.2xlO^^ erg s-^ 



3 ANALYSIS AND RESULTS 

3.1 Constraining Star Formation 

The measurement of the radio angular sizes of HDF 147 and 
HDF130 were made using the AIPS software, JMFIT. JM- 
FIT was used to fit a single 2-d Gaussian to the core com- 
ponent of each galaxy in the CLEANed MERLIN plus VLA 
radio image restored with a 200mas circular beam. Subse- 
quently the restoring beam of this image was deconvolved 
from these fitted sizes. The radio images of both HDF130 
and HDF 147 are of high signal-to- noise and close to the 
pointing centre, thus the measurement of deconvolved ra- 
dio sizes several times smaller than the rest ored C LEAN 
beam of the image can be reliably made (e.g. ICondon et al.l 
Il998l l. The fitted sizes, deconvolved with the image restoring 
beam, of the cores of HDF147 and HDF130 are <80mas and 
<65mas respectively, with an error on the fitted FWHM of 
20mas (Condon 1997). At these redshifts, these angular size 
limits correspond to linear sizes of <640pc and <520pc. 
Figure 1 shows contoured images of HDF 147 and 

ing Facility (STECF/ESA) and the Canadian Astronomy Data 
Centre (CADC/NRC/CSA). 



HDF130 overlaid on ACS/NICMOS imaging. Muxlow et al. 
(2005) classify HDF147 (J123623+621642) as an AGN can- 
didate and note that it has a faint one-sided radio emission 
extending ^-^0.6^^ to the south of the compact core. This core 
and jet-like emission can clearly be seen in Figure 1. The 
radio structure of HDF130 (J123617+621540) is compact 
in the high resolution images presented here, however from 
lower resolution (0.5^^) imaging of these data Muxlow et al. 
(2005) showed this source to have compact core with two- 
sided radio emission extending over a linear extent of 1.3^^ 
The HDF 147 radio morphology of Fig. 1 exhibits a low lu- 
minosity (~20/iJy) extension along the galaxy's minor axis. 
This is a possible candidate for radio jet outfiows but pro- 
vides no obvious clue to source orientation. 

If star formation were the source of radio emission, 
both sources would exceed the maximum star formation 
density possible for a given mass and dynamical timescale 
dElmegreenj Il999l ). HDF147 and HDF130 have limits on 
their radio size smaller than most local ULIRGs. As local 
ULIRGs with T^sfr ~2OOM0 yr~^kpc~^ are already form- 
ing stars near a theoretical maximum (e.g., iTacconi et al.l 
2006), HDF147 and HDF130 would exceed this limit by fac- 
tors of 50-70 X given that their radio- inferred SFR densities 
are 15OOOM0 yr~^kpc~^ and IIOOOMq yr~^kpc~^ respec- 
tively (with only small variations due to geometry and gas 
fractions) . We therefore conclude that the radio emission we 
observe must be largely emitted by AGN. 

We further constrain their star for mation rates usin g 
the rest-UV luminosities (Equation 1 of iKennicuttI 1 1 998l V 
Without removing the effects of dust extinction, we con- 
strain HDF147 to a maximunQ SFR of 5tt^o ^^"^ 
HDF130 to an SFR of 17^^ Mq yr-^ Corrected for ex- 
tinctions of E(B-V)=0.40 and 0.03 respect ivel}|fl (calcu- 
lated from HYPERZ stellar population model fits, as de- 
scribed in the next section), the UV SFR constraints are 
SFRHDFi47=nOtHMQ yr-^ and SFRjfi,Fi30=22t^ Mq 
yr~^. The low 24/im and 850/im luminosities are consistent 
with these SFRs (Ssso < 2mJy; SFR < 3OOM0 yr"^). 



3.2 Characterizing the Stellar Populations 

To learn about the stellar components of HDF 147 and 
HDF130, we fit the optical and Spitzer — IRAC photomet- 
ric points to stellar population te mplates from MarastonI 
(2005) using the hyperz package (|Bolzonella et al.ll2000l V 
HYPERZ output includes internal extinction factor Ay and 
lesi-K band absolute magnitude, which we use to calcu- 
late extinction-corrected UV SFRs and derive stellar masses 
respectively. The Maraston model populations assume a 
Salpeter IMF and range in age from lOMyr to 15 Gyi0- We 
find that the older models (10-15 Gyrs) are much better fits 
to both galaxies' photometry than a 'young' model aged 



^ If all of the UV emission were from star formation. 
^ We considered the effect of IGM extinction at z~2, but found 
it was insignificant compared to magnitude uncertainty and ex- 
tinction internal to the galaxies. 

^ Although the age of the Universe is only ~13.7Gyrs, the aging 
of the Maraston stellar population models has not been normal- 
ized to the current cosmology, thus allowing "15 Gyr" old models. 
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is slightly elevated at z^ 2 than locally (e.g., IPeng et aP 
I2OO6I : iMcLure et alllioOGl ) although such variation could be 
caused by a bia s towards large black holes fcauer et al ]l2007l : 
[Alexander eTal . 2008b). 

Figure [2] explores SED fitting for these galaxies from 
the rest-UV wavelengths to the radio. The SED shortward 
of 1/im is dominated by the best-fit hyperz Maraston stel- 
lar population models (dotted lines in Fig. [2]). Towards 
longer wavelen gths, we int roduce a template Mrk231 SED 
(adapted from lBertall2006l l including both AGN and mod- 
erate star former components and normalize it to the 24/xm 
flux densities (dashed line in Fig. [2])- The inferred star for- 
mation rate from templates with this 24/im normalization is 
^IOMq yr~^. We then raise the radio synchrotron contribu- 
tion to the SED (but we keep the synchrotron slope, a — 0.7 
the same) to flt the galaxies' bright radio luminosities. 

The optical to mid-IR portion of the SED provides a 
good fit to the data: a sum of old and massive stellar popu- 
lations with moderate (^^lOM© yr star formation. The 
submillimetre fiux density limits at do not place any sig- 
nificant additional constraints on the SEDs. In the rest-UV 
[B kV bands) HDF147 is un-detected. HDF130 exhibits 
a moderate blue excess, with associated SFR=22M0 yr~^, 
which may correspond to the starburst detected in its rest- 
UV spectrum. We emphasize that the total luminosity from 
the star forming component is minor compared to the very 
luminous SED-characterized "l.G/xm-bump" stellar popula- 
tion. 



Fig;ure 1. Radio luminosit y (jChaDman et al.ll2003|). scale sizes 
(lAbraham et al.l l2003l . l2004l ). and stellar mass (IConselice et al.l 
I2OO8I , I2OO5I ) against redshift for radio galaxies in several large 
survey fields. HDF130 and HDF147 are highlighted as black dia- 
monds at much larger radio luminosities, scale sizes, and higher 
masses than are typical high-z galaxies. 



only a few Gyrs. Furthermore, we find that the instanta- 
neous starburst models with 0.67 metallicity fit best, where 
goodness of fit is measured by minimizing the statistic. 

We measure extinction factors Ay = 1.3 for HDF147 
and Ay = 0.1 for HDF130, which corresponds to E{B - 
V)iA7 = 0.40 and E{B - V)i3o = 0.03. These derived val- 
ues are used to estimate the rest-UV star formation rates 
described in the previous section. 

We have used the methodology outlined in iBorvs et al. 

(|2005l l to estimate stellar masses from rest frame X-band 
magnitudes, which we measured from HYPERZ to be -26.67 
for HDF147 and -26.81 for HDF130. We assume a mass 
to light ratio of L^/M = 1.9 because of the galaxies' pre- 
dominantly old stellar populations (t>1 Gyr). Both galaxies 
are exceptionally massive, with log(M/M0) = 11.8to3 and 
11.4t°;^ for HDF147 and HDF130 respectively. The con- 
trast to other z^2 galaxies is seen in panel 3 of Figure [1] 
We use the well studied local M^-M^if relation where M-bh 



= /lOOO ( Magor rian et al.l ll998l: iKormendv Richston 3 
I1995I : iGebhardt et al 



2OOOI ) 



8 Q+O-6 
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to estimate black hole masses 
and 8.51q"5, with the caveat 
in the Mbh /M^ 



of log(Msif/M0) 

that there is a factor of two variation 
ratio. We caution that M^-Msij has a larger spread and 



3.3 Size and Morphological Classification 

In contrast to optical imaging, NICMOS near-IR imaging 
reveals that the galaxies have half light radii of 8kpc and 
4kpc, over 20 times the scale of the radio core emission. 
This is shown in Figure El Van Dokkum et al. (2008) find 
that while local galaxies have a typical size of R1/2 '^Skpc, 
elliptical galaxies of ~3xl0^^ Mq at z~2.3 are much more 
compact with average Ri/2=0.9kpc, which is significantly 
smaller than HDF147 and HDF130 (see Figure [T] panel 
2). However the stellar densities of HDF147 and HDF130 
(p.,i47=1.2l J:^ Mq kpc-^ p,,i30=5.7l^:^ Mq kpc'^ 0) are 
much lower than the high-density z^2.3 galaxies, and sim- 
ilar to local ellipticals (Blanton et al. 2005). A similar size 
comparison o f massive _(M >J_O^^Mq) high-z (1.7<z<2.0) 
galaxies from lTruiillo et aL 1 (|2007l ) also shows that HDF147 
and HDF130 are much l arger than typical q uiescent galax- 
ies at similar redshift s. iTruiillo et al.l (|2007l l find that for 
galaxies with large Sersic index (n>2.5) the expected size is 
^i/2=l-li0.4kpc and for shallower potentials (n<2.5) it is 
i?i/2=2.3±1.8kpc. 

We fit Sersic profiles to H-band imaging to search for 
substructure, derive Sersic indices, and model their mor- 
phologies. Fitting both HDF147 and HDF130 using Galfit 
(^Peng et al. I I2OO2I ). we find that HDF147 closely resembles a 
de Vaucouleurs profile (n=3.5±0.3) while HDF130 fits weh 
to a more extended exponential (n=l. 5^0.5). Clearly by this 
analysis, HDF147 is a giant elliptical galaxy with smooth 
morphology and Re '^9kpc. HDF130 is near the boundary 
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Figure 2. A composite SED for the two galaxies (red for HDF147, bl ue for HDF130) i n rest wavelength, from UV through the infrared, 
submillimetre, and radio. The old stellar populations, adapted from Maraston (2005) stellar evolution models, are fit to the optical 
and Spitzer — IRAC photometry, around the rest-1.6/im stellar bump (dotted lines labeled 'M05 Elliptical Models'). A Calzetti dust 
extinction law (e.g. Calzet ti et al.lll99 4') is used to extinct the flux in the UV, with derived reddening factors of E{B — V) = 0.40 and 
E{B -V) = 0.03 for HDF147 and HDF130 respectively. A Mrk231 template SED ( Berta 2006), scaled to the 24/im flux densities, 
is added to represent the far-IR and radio portion of the SED. We flt the SED to the radio points by scaling up the emission from 
synchrotron emission {a = 0.7) from the nominal flux densities that were expected from the Mrk231 template (dashed line). 
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Figure 3. MERLIN radio contours overlayed with ACS B, V, i, z (F435W, F606W, F814W, F850LP) and NICMOS H (F160W) imaging 
of galaxies HDF147 (top) and HDF130 (bottom). HDF147 shows a faint extended region to the north of the radio and NICMOS centers; 
its NICMOS size is substantially larger than most galaxies at its epoch at a ~17kpc diameter. HDF130 has compact blue emission which 
is coincident with the radio center, but has a much more extended morphology, ~8kpc across, in if-band. The ACS panels show only 
the 15/iJy radio contour for clarity while the NICMOS panels show more complete contours: 15/iJy, 40/iJy, 150/iJy, 255/iJy, and 380/iJy. 
The peak flux density of HDF147 is 423/iJy and for HDF130 it is 162/iJy. 



between elliptical galaxies and disk-like structures at n=1.5 
and shows no substructure; its giant size, giant mass, lack of 
significant star formation (<1OOM0 yr~^), and old stellar 
population are characteristic of large quiescent galaxies. 



3.4 Environment Characterization and Further 
Classification 

Given the large stellar masses, old stellar populations, and 
large scale sizes of these galaxies, the question of their 
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environment naturally arises since locally giant elliptical 
(gE) galaxi es reside e x clusiv ely in highly clustered envi- 
ronments. Blai n et alJ (|2004l ) observed that there was a 
strong association of 5 SMGs in this field at the redshift 
z ~ 1.99 (around HDF130), represent ing the largest cluster 
of SMGs in theiChapraan et al.l (|2005l ) spectroscopic survey. 
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(|2008l ) have studied this structure further 
showing that an over-density of galaxies of all types exists 
at this redshift, although not with the same degree of con- 
trast as the SMGs. No such structure was found associated 
with HDF147, although it is possible that one exists and has 
not bee n observed due to s pectroscopic incompleteness. 

The lBauer et al.1 (|2002h study of extended X -ray sources 
detects elongated X-ray emission roughly centered on 
HDF130 (CDF-N source 2, scale size see their Figure 2 
most upper right panel). Due to its double lobed shape and 
sharp cutoff on the ends, the emission could be due to in- 
verse Compton scattering of the CMB due to AGN outflow 
in HDF130 rather than the cooling gas of an extended clus- 
ter. The size of the emission region is 300-500 kpc, and no 
radio lobes or plumes are observed to correlate with the X- 
ray region. Detailed analysis of the X-ray observations of 
this region are presented in Fabian et al. (2009, in prep). 

The se galaxies we r e not selected as z = 2 BX/BM ob- 
jects by ISteidel et al.l (|2004l l from their colo rs, and they 
do no t appear in the spectroscopic catalogs of ICowie et al.l 
Both galaxies satisfy the p-BzK selection criteria for 
passively evo lving galaxie s (BzK < -0.2, z — K > 2.5) as 
described bv iDaddi et al ] (|2004l ). However, the scale size of 
these galaxies is much larger than the vast majority of p- 
BzK galaxies which are generally very compact (scale radii 
<0.9 kpc; .van Dokkum et al...2008. ). 



jjjy RADIO GALAXIES AT z~2 



SUBMM-FAINT 
STAR FORMING 
RADIO GALAXIES 



SUBMM 
GALAXIES 



SFRGs 30% 



SMGs 25 

SMGs WITH AGhT 

8% 




3% 

T 

HDF147. 
HDF130 

EVOLVED GALAXIES 



Figure 4. A schematic pie chart characterizing the z~2 /iJy 
radio source population. SMGs and SFRGs have optical spectra 
consistent with starbursts, while AGN and SMGs with AGN have 
spectral signatures of AGN. The SMGs and SMGs with AGN 
are detected at 850/im. HDF147 and HDF130 are highlighted as 
the small red slice, between SFRGs and AGN; while they were 
previously diagnosed as SFRGs from their spectra, we have shown 
that they are more appropriately classified as AGN. 



presence of an AGN, low X-Ray luminosities, no evidence 
for high ionization UV emission lines and low 24/im lumi- 
nosities. 



4 DISCUSSION 

To understand the importance of galaxies like HDF147 and 
HDF130 at z^2, we first characterize the parent sample 
of radio sources from which they were selected. We illus- 
trate schematically in Figure |4] the relative fractions of dif- 
ferent galaxies (SMGs, SFRGs and AGN) comprising the 
/iJy z~2 radio galaxy population (a more detailed analy- 
sis is in Casey et al. 2008, in preparation). HDF147 and 
HDF130 are highlighted since our observations have demon- 
strated that they should be classified as AGN-dominated 
rather than star formation-dominated SFRGs. These two 
galaxies represent a rare type of massive galaxy with an esti- 
mated volume density of 1.5x 10~^ Mpc~^ (given two objects 
detected in a well studied area 15^x15^ and redshift range 
1.8<z<2.6). Evolved elliptical galaxies of this type (with 
stellar masses over 10^^ Mq and ne^li^ibl e black hole ac- 
cretio n) are extremely rare at high redshift. IConselice et al.l 
(|2007l ) estimate $ < 5x10"^ hfo Mpc"^ the volume density 
estimate for highly evolved z^2 gal axies, which is less than 
$ for SMGs (jChapman et al.1 120031 ). While this suggests a 
quite sparse population, a combination of spectroscopic in- 
completeness in the radio source population and a possible 
wider range of physical properties consistent with these two 
galaxies could suggest that they are an important evolution- 
ary stage of massive galaxies. 

This section explores physical explanations for the ob- 
servational constraints: high-resolution radio indicating the 



4.1 Evolved Galaxies with Poorly Accreting Black 
Holes? 

The observational constraints (low 24/im and X-Ray lumi- 
nosities) require that there be very little accretion (77 ^0.1) 
onto the massive 10^-^ and 10^-^ Mq black holes. This lower 
accretion efficiency could be explained by decreased fuel- 
ing of the AGN, where most cool gas has been converted 
into stars or expelled from the elliptical galaxies by feedback 
mechanisms. A lack of hot AGN torus dust is a condition for 
this interpretation because of the weak mid-IR luminosity. 

With low X-Ray luminosities yet high radio luminosi- 
ties, we suggest these radio fluxes may be boosted by some 
mechanism. We propose weak radio beaming, here deflned 
as jet- line of sight orientation angles less than 30° with 
fluxes boosted by a factor of a few, to explain the 200- 
500/iJy radio flux densities. At the other extreme, refer- 
ring back to Figure (2] we see that the radio emission in- 
ferred from a moderate-luminosity AGN with modest star 
formation would be ^2/iJy. The intrinsic emission of the 
HDF147/HDF130 AGN components must lie between 2 and 
200/iJy. 

We suggest that HDF147 and HDF130 are likely 
the high- redshift ana log of local Fanaroff- Riley Fs 
( iFanaroflP RilevI I1974I I. which we use to estimate 
volume density and the probability of our observations. 
Both galaxies fall below the critical FRI/FRII boundary 



Two z^2 MERLIN Evolved Radio Galaxies 7 



in radio luminosity and rest-frame optical magnitude 
(i.e. the r adio-optical magnitud e bivariate FRI/FRII 
boundary; Ledlow OwenI Il996l ). They are also classi- 
fied as FRI's in the original classification scheme (i.e. 
L < LcRiT (178 MHz) ^ 2x 10^^ WH z~^ str~^ assum- 
ing a = 0. 8; iFanaroff RilevI Il974h . Their intrinsic 
emission is likely similar to the average luminosity 
of these nearby FRI AGN in evolved galaxies with 
uKilAGHz) -5x10^^ erg s"^ The probability of observ- 
ing intrinsic radio luminosities of vLu (1.4 GHz) ^ 10^^'^ erg 
s~^ a mong FRI radio sources is very low ([Ledlow OwenI 
Il996l ): 0.5% (to observe HDF147 at = lO^^"^) and 1.3% 
(to observe HDF130 at z/L^ = lO'^^-^). Both probabilities 
represent the fraction of FRIs in the Ledlow & Owen 
sample with that luminosity or greater. If we allow beaming 
in both cases, so that intrinsic luminosities drop by a factor 
of ten, then the probability of our observation increases 
marginally, by 5 times. We cau tion that the statistics 
behind the iLedlow k OwenI (|l996l ) sample are limited and 
only include nearby sources, and that the evolution with 
redshift might not be constant. 

We can also use the much larger and well studied statis- 
tics of the faint-end quasar luminosity function, and apply 
constraints based on different beaming-strength geometries. 
The most recent studies of the quasar lumi nosity function 
(QLF ) are from iHopkins et al.l (|2007l l and 'Ric hards et al.l 
((20061), whose faint-end reliabilities, as they relate to proba- 
ble low- luminosity low- accretion rate AGN, are discussed in 
ICasev et all (|2008l ). For 26.5<i<23.5 (the magnitude selec- 
tion cut on the parent sample of SFRGs), and 1.7<z<2.3, 
the AGN surface density is 168=b2 deg~^. Assuming an 
isotropic distribution in orientation angles for type 1 AGN, 
the number of these faint AGN whose jets are pointed to- 
wards Earth within 30^ (weak-beaming) is 42±17 deg~^. 
For strong beaming (<3^ opening angle, ^-^100 amplifica- 
tion) the numbers shrink to 0.8±0.3 deg~^. We therefore 
expect 1.2±0.5 low- luminosity, weakly beamed AGN to ap- 
pear in the MERLIN coverage of GOODS-N, which agrees 
with the two galaxies we have found. 

Since we classify these galaxies as massive, with old stel- 
lar populations, there is also a constraint on X-Ray emis- 
sion produced by hot gas that we test. While all massive 
ellipticals in the local Universe have large quantities of X- 
Ray emitting gas ( as have been shown out to 2; = 0.7 in 
Lehmer et al. I2OO7I ). the expected X-Ray luminosity from 
hot gas in these systems is about an order of magnitude be- 
low the current Chandra X-Ray limits (Lx'^10^^ ergs~^ jf|. 

4.2 Alternate Hypotheses 

The possibility exists that HDF147/HDF130 are highly- 
beamed (opening angle <3^) low- luminosity AGN. This 
hypothesis would further reduce the concern for low dust 
masses in a hot torus since low- luminosity AGN would not 
create significant re-radiated mid-IR fiux even in the pres- 
ence of dense dust clouds. However, looking straight down 

^ The hot-gas X-Ray constraints were calculated using the av- 
erage X-Ray to K-band luminosit y ratio of early type galax- 
ies from El lis O'SullivanI (|2006l V iLehmer et"al1 (,20071 and 
iLehmer et al.l (l2008l V 



the radio jets, one would expect signs of a BL Lac type 
optical spectrum, which is clearly not present. A variation 
on the low- luminosity AGN hypothesis is that the central 
black holes are much less massive than the local M^-Mbh 
predicts. The lower mass limits determined by Eddington 
limited accretion are 10^'^ and <10^"^ Mq (two orders of 
magni tude below the predictions of M^-Mbh, as in Figure 
3 of Borv s et al.1 12005). However, under this hypothesis we 
would have to advocate factors >100x beaming to boost 
the radio lumino sity above what is e xpected for such small 
black holes (e.g.. iMerloni et al.l l2003'). With the analysis of 
the previous subsection, the orientation geometry is likely 
implausible for either variation of the highly-beamed, low- 
luminosity AGN hypothesis (with probability ^1%, requir- 
ing both galaxies to have jet-line of sight orientation angles 
of <3"). 

Compton thick AGN are also an alternate hypothesis. 
The lack of detection in the soft X-ray could be explained 
by a deeply buried, obscured AGN. The extreme obscuration 
could also explain the weak rest-UV. The AGN would per- 
meate the dust obscuration in the radio, however one would 
expect a much higher 24/im fiux density d ue to a high dust 
content (by at least an order of magnitude I Alexander et all 
l2008bl ). 



5 CONCLUSIONS 

We have presented mult i wavelength observations of two high 
redshift galaxies, arguing that they are highly evolved giant 
elliptical galaxies with low- luminosity, beamed AGN. Their 
near-IR morphologies (NICMOS H-band) show giant ellipti- 
cal characteristics typical of an old stellar population whose 
stellar bump is well characterized in the IRAC filters. The 
physical sizes of both galaxies are 2-5 times larger than most 
evolved galaxies at the same redshift, and similarly the stel- 
lar masses are large, around '^4x10^^ ^0- M^-Mbh 
relation, we infer black hole masses '-^3xlO^M0. The lack 
of strong AGN suggests accretion at very low rates. Without 
X-Ray detections and without strong power law continua in 
the rest-UV, we assume they are low- luminosity AGN whose 
radio emission must be beamed by a factor of a few. We offer 
discussions of other possible scenarios, but it is most likely 
that we are observing evolved galaxies with AGN that have 
nearly exhausted their fuel source. While this type of highly 
evolved galaxy is rare at high-z, the supplemental knowledge 
that galaxy over-densities have been found near z = 1.99 in 
the HDF reinforces, if only slightly, the notion that galaxy 
evolution occurs more rapidly in cluster environments. 

These galaxies were originally classified as radio- 
selected star formers, likely with hotter dust temperatures 
than their SMG counterparts because of their 850/im non- 
detections. Only their compact MERLIN radio morphologies 
and weak 24/im fiux set them apart from high-SFR galax- 
ies. A full analysis of the SFRG population is in Casey et al. 
(2008), in preparation. The planned eMERLIN deep point- 
ings at 1.4GHz and 5GHz in GOODS-N wih push detections 
to '^0.5/iJy rms, measuring radio spectral indices and spa- 
tially separating emission from AGN and star formation. 

So far these galaxies have gone unnoticed in an other- 
wise well studied deep survey field. The potential to study 
more galaxies exhibiting the same properties is possible. 
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as many have probably gone un-noted in studies due to 
incredibly faint X-ray and optical fluxes. Although both 
objects were measured at 1.4 GHz with the VLA, it was 
the MERLIN high resolution coverage which distinguished 
the extreme radio luminosities as emanating from the ex- 
tremely compact galaxy cores. If HDF147 and HDF130 con- 
tain beamed low-luminosity AGN as we propose, they pro- 
vide a testbed for understanding evolved and poorly accret- 
ing systems at high redshift. 
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